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A RESISTOR-NETWORK SOLUTION FOR THE ELASTO- PLASTIC 
TORSION PROBLEM 


J. H. Weiner,! M. G. Salvadori,2 M. ASCE, and V. Paschkis? 


INTRODUCTION 


The solution of torsion problems in which the stresses exceed the elastic 
limit presents great analytical difficulties. Numerical solutions, using relax- 
ation techniques, have been developed? and are effective. However, if many 
mesh points and separate calculations are required, as for example in obtain- 
ing an accurate torque vs. angle of twist relationship, the relaxation technique 
becomes quite laborious. 

In this paper a method of solution using a resistor network is presented 
for the torsion of a cylinder of perfectly elasto-plastic material with convex, 
simply-connected cross section. The finite-difference approximations cus- 
tomarily used in the relaxation procedure are made, but the relaxation is 
accomplished automatically by the network. Thus, once the network for a 
given geometry is set up, a complete torque vs. twist curve may be obtained 
in a short time. 


Electrical Analog of the Elasto-plastic Torsion Function 


The method will be described as applied to the torsion of a square section 
of side a. The stress-strain relationship of the material is assumed to be as 


shown in Fig. 1. 
In the usual notation, the components of shear stress on the cross section 
are given 
=- o¢ (1) 


Tax = $F; Tay ~~ Ox 


where the stress function ¢(x,y) satisfies the equation: 


A244 2G6 =0 (2) 

in the elastic portion of the cylinder and the over-riding conditions 
|44|= Toi (3a) 
+ 2G0 90 (3b) 


in the plastic portion. Equation (3a) corresponds, in terms of the membrane 
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analogy, to the erection of a roof of slope To confining the membrane. Equa- 
tion (3b) corresponds to the restriction that the roof may only push, not pull, 
on the membrane. In addition ¢ satisfies on the boundary C the condition: 


%, = 0 (4) 
The torque, M, corresponding to a given $(x,y), is given by: 


M = 2S S dx dy (5) 


For the study of elasto-plastic torsion it is convenient to introduce the 
following dimensionless quantities: 


§ = ; N= 
Then equations (1) to (5) become, 
(1') 
(2') 
(3'a) 


(3'b) 


M = 2Ga° Swat dn = 2G° (5") 


Introducing the usual finite difference approximations, with mesh size 


h = + , Fig. 2, equation (2') becomes: 
- 4y. 
h2 


In finite difference terms, equation (3'a) becomes, using equation (4'): 


mhrTo 
G m=1,2,....,0 (3"'a) 


where Wn is the value of Wat any node a normal distance mh from the 
boundary. Similarly equation (3'b) becomes: 
+H, - 40, 
h2 
for any node i in the plastic region. 
Consider now the resistor network shown in Fig. 3. Applying Kirchhoff’s 
law to the currents flowing into a typical node, the following relation is ob- 


A“w+ 20a > 0 
Vo = 0 
| 
tained: 
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If E is made very large relative to any e., the last term is approximately 
equal to E/Ry and the above equation becomes: 
e +@,+e +e - 4e, 

r = b 1 E 0 as 


If a scale factor between e and W is established 
e = KW, (7) 
equations (2"') and (2"'') become identical for 
(8) 

2R,9 ah? 

The analogs of equations (3", then become: 


(a) e,, Kmh T,/G (3"' 'a) 


e +@€,+e +e - 4e, 


where e_ is the value of e at any node a normal distance mh from the boun- 
dary. 

The boundary condition onW, equation (4'), is satisfied by shorting the 
boundary nodes together and connecting them to ground, i.e., the zero of the 
voltage scale. 

In order to satisfy equation (3"'' a), a low internal impedance voltage 
source is required with taps at mKh T,/G (m =1,2,...,n). In practice, K 
is chosen so that Kh T,/G = @g, the voltage across a single cell of a storage 
battery. Starting at the node closest to the centroid of the section, a normal 
distance h from the boundary, since this is the node with highest slope, the 
voltage is compared with e, by means of a difference amplifier. If the nodal 
voltage is greater than es, it is shorted to e,; if it is lower than e, it remains 
untouched. Then the next nodes in the same row are compared with e, in an 
identical way. Following examination of the nodes in the first row, the nodes 
in the second row a normal distance 2h from the boundary are compared with 
2e., etc. 

It is seen that equation (3''' b) is automatically satisfied for any node com- 
pletely surrounded by nodes in the plastic region. A numerical check, made 
by substituting the appropriate voltages in equation (3"''b) is required for the 
nodes, on the border between the elastic and the plastic regions, and those at 
which it is not satistied are disconnected from the battery. However, if the 
nodes are connected to the battery in the order described above, this reversal 
will seldom occur. 

In this manner, there is established a complete analogy between the finite 
difference equations governing W and those governing the nodal voltages in 
the electrical circuit. To determine the values of W, therefore, it is only 
necessary to measure the nodal voltages and use equation (7) to get the equiv- 
alent values of W. 

To determine the torque, M, corresponding to these nodal values of Y, it 
is convenient to use the trapezoidal rule for numerical integration of V equa- 
tion (5'): 2 

(2n-1) 3,2 (2n-1) 
M = 2Ga3v 2Ga9n? 5") 
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where the summation is carried out over all the internal nodal values. It 
should be noted that it is theoretically possible to determine the torque with a 
Single current-difference measurement. For the current I flowing into the 
circuit from the voltage source E for a given set of voltages e; is: 


(2n-1)2 (2n-1)? 


where I. is the current which would flow into the circuit if all the e; were 


zero. Therefore: 


2G an‘r,? 
K 


(lo - = 


However, since I, - I is very small relative to I, its measurement is deli- 
cate and further work is required for the development of a technique giving 
sufficiently accurate results. The values of M inthis paper, were obtained 
by means of equation (5"). 


M = (I, - 1) Goat (5""") 


Determination of Torque vs. Twist Curves 


It is most convenient to present a complete torque vs. twixt curve in terms 
of the dimensionless quantities M/M, and 6/@, where M, and 9, are, 
respectively, the torque and angle of twist at which plastic deformation just 
begins, i.e., at which the maximum stress in the cross-section is T.. 

The dimensionless slope of the roof is therefore first determined from an 
“elastic” network experiment, that is one in which the over-riding condition, 
equation (3"), is not applied. For this experiment, any convenient value of K 
may be used. The maximum value of vw , occurring at the boundary point 
nearest the centroid of the section, may be computed from the nodal values of 
W by the use of suitable finite-difference formulae. 

The above experiment yields the values of M,, 6, and the value of T/G = 
| Vw) max. The value of K is now chosen so that Kh %/G = @g, the voltage 
across the single cell of a storage battery, and this choice of _K is seen from 
equation (8) to be equivalent to a choice of E, corresponding to the above 
value of 4,. From equation (8) it is further seen that if the scale factor K is 
to be maintained constant, the imposition of an angle of twist 6 = m 6,(m>1), 
corresponds to increasing E to mE>. 

The technique described above is then applied for E = mE, to obtain the 
torque M corresponding to an angle of twist m6,. 


Experimental Results 


The technique explained in the previous sections was used to determine 
the torque vs. twist curve for a cylinder of square cross-section on the Heat 
and Mass Flow Analyzer at Columbia University School of marries. This 
equipment is suitable for circuits with as many as 576 nodes. »7 The three 
curves of Fig. 4 give M/M, vs. 6/6, for h = 1/6 and h = 1/12 with a roof 


6. V. Paschkis and H. D. Baker, “A Method for Determining Unsteady State 
Heat Transfer,” Trans. ASME, Vol. 67, (1942) p. 105. 

7. V. Paschkis and T. Mirsepassi, “Temperature Distribution in the Hearths 
of Blast Furnaces,” Iron and Steel Engineer, 1954, pp. 53-66. 
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slope A approximated from the nodal voltages of the “elastic” network experi- 
ment by the first difference expression at the boundary: 


(4 


Kh ’ 


and for h = 1/12 with the roof slope approximated by: 


A" = (Ae - sare + 74 3e---), (9) 
The differences higher than the second of the “elastic” voltages e become 
erratic in the experiments and the best slope for the accuracy used in these 
measurements is obtained by the first two terms of equation (9). It will be 
noticed that the results are very sensitive to the value of A, while they are 
not substantially influenced by lumping. 


Error Evaluation 


In order to check the accuracy of the experimental results some of these 
were checked against results obtained by relaxation. 

Table I gives the experimental and relaxation values of W4,a for 9/6, 

1.5 and h = 1/6 with the slope A'. The largest and average errors equal 1. Dh 
and 0.87%, respectively, and are to be compared with the 1% error in the 
voltage measurements. It is thus assumed that the other experimental re- 
sults are also affected by errors of the order of 1%. 

The value of the double integral V in equation (5") for 46 /9, = 1.5 com- 
puted by the trapezoidal rule from the experimental values has an error of 
only 0.05% relative to the relaxation values. The same error in the value of 
the integral for the elastic case equals 0.33). The experimental value of the 
ratio M/M,, which depends on the elastic case integral, has an error of 
- 0.30% measured against the relaxation results. 

The purely numerical techniques used to evaluate the pivotal values of W 
and M introduce three sources of errors inthe results: a) the approximate 
evaluation of the pivotal values of W due to lumping, i. e., due to the spacing 
h; b) the approximate evaluation of V , equation (5"'), i. e., the volume under 
the analog membrane; c) the approximate evaluation of the boundary slope A . 

To evaluate the influence of lumping it may be noticed that the error in M 
is of the so-called h@-type and an h*-extrapolation may be used from the 
values for h = 1/6 andh = 1/12. For the elastic case, the extrapolation gives: 
My = 0.1377 - 0.1291] G6, at = 0.1406 G6,a4, 
which checks up to the last figure with the result given in Timoshenko® and 
indicates an error of 8.18% in the h = 1/6 approximation and of 2.06% in the 
h = 1/12 approximation. For the elasto-plastic case 6/6 = 1.5 the extrapo- 
lation would not be valid since the h = 1/6 andh = 1/12 cases are evaluated 
by means of different roof slopes, but it may be presumed that the same or- 
der of magnitude of error as in the elastic case will be present here, i.e., 
that M be in error by approximately 2%. 

In order to evaluate the influence of the integration procedure, V in equa- 
tion (5"') was computed by the trapezoidal rule and by Simpson’s 1/3 rule on 


8. S. Timoshenko, “Theory of Elasticity,” McGraw Hill Book Co., New York, 
1934, p. 249. 
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the basis of the relaxation results (h = 1/6) for the elastic case and for the 
case 6/6, = 1.5 with the slope A'. It was found that the difference between 
the two evaluations of the integral equals 6.20% in the elastic case, and 4.98% 
in the case 6/6, = 1.5. On the other hand, the ratios of the integrals V 
(which is identical with the ratio M/M,) evaluated by the two methods, has a 
difference of only 1.3% for 6/6, = 1.5 when the slope is A' and of 0.71% when 
the slope is A". 

Since the error in Simpson’s rule is of order n4 it may be concluded that 
the error in the ratio M/M, due to integration by the trapezoidal rule, which 
is of order h”, for h = 1/6 is less than 1%. 

To compute the error due to the roof slope, the integral V of equation (5") 
was evaluated from the relaxation results for h = 1/6 and 6/6, = 1.5 with the 
two values A' and A" of the slope, using the trapezoidal rule. The error in 
the A’ solution, assuming the A" solution to be correct, equals 17.18% and is 
identical with the error in M/M,- The same error using Simpson’s rule 
equals 17.66%. 

It is thus seen that the errors due to other causes are negligible compared 
to the errors due to the approximation of the roof slope and that this particu- 
lar parameter should be evaluated as accurately as possible whatever the 
method of solution adopted. 
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Node 


| Network | Computed 


| 
1776 1793 


1433 
1583 1585 
641 
865 
865 


Nodal values of a 
1/6, Slope A', = 1.5 
See Figure 2 for location of nodes. 


a 
Table I 
| 
| 
% Error 
0.95 
| 2 | +0.35 | 
3 | -0.13 | 
4 +1.1 | 
5 
6 +1.2 
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Fig. 2 Finite Difference Mesh 
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Fig. 3 Resistor Network 
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